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Although a discovery of wimps either at colliders or in direct experiments would have enormous 
implications for our understanding of particle physics, it would imply less than one would like 
about our understanding of the dark matter in the universe or in the galactic halo: it surely is 
possible that the discovered particles account for only a little of the total dark matter. To establish 
the cosmological significance of a wimp discovery, their density must be determined. We show 
that data from neither hadron colliders nor direct detection experiments alone can be sufficient 
to determine the local or relic density of discovered wimps, even allowing all needed assumptions 
about cosmology and astrophysics. We provide examples of dark matter models where nearly 
identical detector or collider signals correspond to very different densities. We show, however, that 
it may be possible to determine the density of wimps by combining data from both experiments; 
we present a general method to do this in the case of supersymmetric dark matter, and describe 
how similar studies could be made for other wimp candidates. 
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I. INTRODUCTION 

We are now confident that our universe contains a large amount of cold dark matter. The most popular 
particle candidates for dark matter are stable, weakly interacting massive particles (wimps). These particles 
are being searched for directly and indirectly by many experimental groups throughout the world, and may 
be produced at colliders in the near future. But we must ask, what is the cosmological significance of a 
discovery of wimps? The answer depends on what fraction of the total dark matter they compose. 

If wimps are discovered experimentally, there are many reasons to hope (and perhaps expect) that they 
solve the dark matter problem — i. e. that the relic density 1 of the wimps is equal to the known relic density of 
cold dark matter. A universe in which almost all the dark matter was composed of a single kind of particle 
would be one arguably kind to physicists. But of course, whether or not we live in such a universe is a 
question that can only be asked and answered experimentally. 

Throughout our discussion, we will find that the connection between the density of discovered wimps and 
the data seen at colliders or in direct detection experiments is quite tenuous: it will be all too easy to find 
models which agree with the initial data having almost any desired relic density. To only consider models 
consistent with the data that solve the dark matter problem is to avoid the physical question of the nature 
of dark matter. To use the known relic density of dark matter to constrain the relic density of a particular 
wimp candidate is to make unjustified claims based on wishful thinking. Put another way, cosmological data 
does not imply a lower bound on the relic density of any particular wimp; there exists only an upper bound: 
no wimp can compose more than the entirety of dark matter. 

The question of the composition of dark matter is usually avoided by the implicit assumption that only 
one wimp would exist in nature. Of course this hope would need to be verified experimentally. But even if 
there were only one wimp particle in nature, there would still be no reason to assume that its relic density 
were equal to that of cold dark matter. This is because wimps are not the only kind of particle that can 
contribute to the total dark matter of the universe. For example, there are strong theoretical motivations 
for the existence of axions; if axions exist, they may account for a significant fraction of the total cold dark 
matter. But it may not be possible to accurately compute the amount of axions produced (non-thermally) 
in the early universe 0. It may be that the only way to determine what (possibly significant) fraction of 
dark matter is axionic is to determine the relic abundance of all other contributions. Therefore, the only way 
to solve the dark matter problem, will be to use experimental data to 'measure' the actual relic abundance 
of discovered wimps. As we will find, this could only be (and perhaps can be) done using data from both 
colliders and direct experiments combined. 

Although the production of wimps at colliders would have enormous implications for our understanding 
of elementary particle physics, it implies less than one would like about dark matter. As it turns out, it may 
be exceedingly difficult to calculate their (thermally-produced) relic density with data from hadron colliders 
for general wimp models. This is even true for specialized models, such as supersymmetric dark matter. 
We give examples of models with nearly indistinguishable collider signatures that have vastly different relic 
densities. Furthermore, there is no reason to expect that all the dark matter was produced thermally, or 
that thermal production took place 'canonically' in the early universe — e.g. in the absence of quintessence. 
We must therefore directly measure their density in the galactic halo, and this cannot be done at colliders. 
But even if these issues could be dealt with, the needed particle physics information is unlikely to exist. 



1 Throughout this paper we will try to distinguish between the cosmological and local (halo) relic densities of wimps. Although 
both of these are 'relics' of the early universe, we will exclusively use 'relic density' for the 'cosmological' relic density. In 
Section II.C.l, we describe why such care is necessary. 
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Fortunately, data from direct detection experiments is directly related to the actual local density of wimps 
in the galactic halo. However, to extract the local density from these experiments, the wimp must be iden- 
tified and its couplings to matter must be calculated. And this information is only obtainable at colliders. 
Until the wimp is observed and studied at colliders, the local density will be invisible to direct experiments. 
These issues, for colliders and direct detection were first raised technically although not resolved in Ref []]. 
Again, we give examples of models with nearly indistinguishable in/direct detection signals that have very 
different relic densities. 

To determine the local halo density of a new particle, its scattering rate with matter must be observed 
in direct detection experiments and its cross sections computed using data from colliders. But this requires 
that the wimp produced in colliders is in fact the same as that observed directly. If the particles observed 
in each had the same mass, it would (strongly) suggest that they were the same. We show that it may be 
possible to determine the wimp mass using data from direct detection experiments alone, and discuss the 
difficulties of this measurement at hadron colliders. 

Once the local density of a new particle is determined, we can deduce what fraction of the local dark mat- 
ter it represents. This will experimentally establish the relevance of the particle to cosmology. Furthermore, 
analysis of this type could lead to novel and strong limits on high-scale physics including the existence of 
additional, high-scale gauge symmetries or quintessence dark energy. 

This paper is organized as follows. First, we study how hadron colliders and dark matter experiments 
can study dark matter independently. We show by explicit examples that, in contrast to what is suggested 
by some specialized studies, it will be extremely difficult (impossible?) to determine the relic abundance 
of generic dark matter wimps produced thermally in the early universe using hadron colliders alone. We 
emphasize that even if the thermal production calculation were possible, there would remain very serious 
uncertainties regarding non-thermal (and non-standard) production. 

We study how dark matter wimps in the galactic halo could be observed directly or indirectly in experi- 
ments. We explain why — and demonstrate with explicit examples — that no combination of data from these 
experiments alone can determine the local density of wimps unless the particle has been identified and its 
scattering cross sections known. To calculate the coupling of a discovered wimp to the matter in detectors, 
one must know many details of its particle physics; of course this requires that it be identified. This infor- 
mation can only be obtained at colliders. We present a general method to determine the local density of 
wimps: using direct detection experiments to measure the product of the scattering cross section and their 
local density, and colliders to identify the particle and estimate its cross sections. We discuss important 
astrophysical caveats of this line of reasoning, and describe possible ways for these issues to be resolved. 

In order to explicitly show how direct detection data depends on the particle physics and local density of 
a wimp, we briefly review how dark matter particles are detected in these experiments. This will allow us to 
introduce formalism common in the literature, and show how data from experiments with different detector 
materials and at different recoil energies can be combined to improve our understanding of dark matter in 
the halo. This discussion will naturally motivate a general (and possibly robust) technique to determine the 
mass of particles seen in these experiments — which is critical to its identification. 

We then consider (in detail) the case where the wimp seen in experiments is the lightest supersymmetric 
particle (LSP) predicted by minimally supersymmetric extensions of the standard model (MSSMs). We 
derive a robust upper bound for (what amounts to) the LSP-nucleon scattering cross section within the 
framework of the most general MSSM using only sparse data from colliders. This analysis applies to the 
most general softly-broken MSSM with a neutralino LSP; no presumptions are made about supersymmetry 
breaking, the masses of sparticles, or mixing phases. We show how this bound can be combined with direct 
detection data to determine an absolute lower bound on the local density of neutralino dark matter 0, Q| . 
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II. THE RELIC AND LOCAL DENSITIES OF WIMPS: COLLIDERS AND DETECTORS 

There is an analogy to be drawn relating our understanding of the dark matter problem to our attempts 
to solve it: the two principal sources of evidence for the existence of dark matter in the universe being 
cosmology and astrophysics; the two ways in which the problem is hoped to be solved being collider physics 
and in/direct detection experiments. 

Studies of the cosmic microwave background and the formation of large scale galactic structure indicate 
that the universe contains an enormous amount of cold (hence massive) matter particles — an amount which 
greatly exceeds that of visible matter in the universe. The excess is called 'cold dark matter' (cdm). Given 
in terms of its ratio to the critical energy density of the universe multiplied by Hubble's constant h, the 'relic 
density' of cold dark matter is confined to the range 0.094 < Vt c ^ m h 2 < 0.129 4]. This is the (universal) relic 
density of all cold dark matter, whatever it is composed of; popular components of dark matter include old 
stellar material, axions, and the lightest supersymmetric particle — all of these and perhaps more making up 
the total density. 

Independent of (and perhaps more forceful than) these cosmological studies, it has been known for more 
than seventy years that vast 'halos' of invisible matter are necessary to explain the anomalous dynamics of 
galactic clusters and galaxies — including our own. For example, studies of the rotation of the Milky Way 
in the vicinity of the sun show that the local density of our galactic halo is roughly ~ 0.3 GeV/cm 3 0, |g- 
This matter in the halo, (perhaps) appropriately also named dark matter, is not confined to the planes of 
galaxies and so must interact at most weakly (and gravitationally) with ordinary matter and itself. This fact 
is strengthened by the non-existence of exotic heavy isotopes, which implies that the halo dark matter does 
not interact strongly or electromagnetically 2 Q. Because both cosmological and astrophysical data imply 
the existence of cold dark matter, the two are supposed to be the same. In short, we know the local density 
of dark matter and we know the cosmological density of dark matter; but the sources of these two facts are 
quite distinct. 

Similarly, there exist two general approaches to solve the dark matter mystery: producing dark matter 
particles at colliders and directly observing them in the galactic halo. Surely it would be very exciting to 
discover wimps in the halo or produce them at colliders. But how does the existence of wimps account 
for the dark matter in the universe? Certainly the mere existence of such particles implies little about 
cosmology. To authoritatively solve the dark matter problem, one must authoritatively demonstrate that 
discovered particles have the correct local or relic densities. It is not sufficient to show that there exists 
a model in agreement with data in which the wimp composes all of the dark matter in the universe; such 
analyses — likely to follow the discovery of wimps — merely avoid the mandate of determining the density of 
discovered particles. Like the case of cosmology and astrophysics, colliders and direct searches offer distinct, 
crucial clues required to solve the puzzle. 



2 If electrically or strongly charged dark matter existed in the halo, it would likely bind to nuclei causing the appearance of 
'exotic' heavy isotopes. 
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A. The (Thermal) Production of Dark Matter in the Early Universe and at Colliders 

The popularity of the supposition that dark matter is composed of wimps is the result of a complicated 
mixture of circumstantial evidence (together with a bit of intuition, hope, and the impressive ease of proposing 
new wimp models — as is evident by a glance at the literature). Experimentally, data from cosmology 
(including large-scale structure formation simulations) indicate that dark matter is massive (cold) and non- 
baryonic; data from astrophysics (including searches for exotic heavy isotopes) indicate that dark matter is 
at most weakly interacting. That dark matter would be well explained by 'wimps' is readily encouraged by 
a rough estimate of what their relic density today would be from thermal production in the early universe. 

In the hot early universe, the production and annihilation of wimps would have been in equilibrium. When 
the universe cooled to an ambient temperature too low for their production, equilibrium would be lost and 
their density would decay exponentially by annihilation. This would continue until the universe expanded 
to such an extent that wimps would no longer effectively 'find each other' to annihilate, and the (comoving) 
density would 'freeze-out.' A back-of-the-envelope calculation for a generic wimp — using a weak-scale mass 
and annihilation cross section and a slightly non-relativistic average velocity — is extremely encouraging: one 
obtains f2 w im P ^ 2 ~ 0.1, the cosmological result 3 . This calculation should not be considered as encouraging as 
it may at first appear — and a negative result should not have been discouraging — because proper calculations 
for explicit models can give relic densities that vary over several orders of magnitude. 

In general, the amount of thermally-produced dark matter in the universe is calculable for any explicit 
model. If colliders observe wimps, then it is this way in which it is (naively) hoped to establish what fraction 
of dark matter they compose: if one knows (or guesses) enough of the particle physics describing the wimp — 
its mass and many of other lagrangian parameters — then one can straightforwardly compute the thermal 
contribution to its relic abundance. 

There are at least two important shortcomings of using hadron collider data to estimate the thermally- 
produced contribution to the relic abundance of wimps: one is practical (and perhaps unsurmountable); the 
other is one of principal. We describe each below in turn. 

1. Why the Thermal Relic Density Cannot be Determined with Hadron Colliders Alone 

Unfortunately, the number of parameters needed to reliably calculate the thermal contribution to the relic 
abundance of wimps is typically so large that there is little hope of ever knowing the result. It should be 
emphasized that in general the relevant parameters are meaningful, physical quantities that are in princi- 
ple measurable once the relevant new physics can be studied experimentally. But the inherent difficulty of 
obtaining the relevant high-precision measurements at hadron colliders makes the task appear unsurmount- 
able: the calculation crucially relies on precise knowledge of too many lagrangian parameters. For example, 
in the case of a supersymmetric wimp (the lightest supersymmetric particle), virtually none of the tens of 
required parameters are known to be measurable in a model-independent way at hadron colliders. Therefore, 
although one would hope to use collider data to determine at least the thermally-produced contribution to 
the relic density, this may be a practical impossibility. 

To demonstrate the loss of generality typically made for the sake of progress, consider the case of 
supersymmetry — certainly the most actively studied model for dark matter. The most general softly-broken 



3 For a slightly more explicit demonstration of this derivation see e.g. Ref. 0. 
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FIG. 1: The relic densities of more than 2500 constrained MSSMs with a Bino-like LSP as a function of m x . The 
mass of the lightest squark is indicated by color. Notice that there is almost no correlation whatsoever between mass 
and relic density. The lack of correlation is clear even when one knows the mass of the lightest squark. 

MSSM lagrangian contains more than one hundred currently unmeasured parameters, many of which enter 
the calculation of the relic density of the LSP. All of these are eventually measurable when superpartners 
and their interactions can be studied experimentally — but a linear collider may be necessary to achieve this. 
To make the problem tractable, most researchers chose to greatly reduce the number of unknown parameters 
by making a number of somewhat ad hoc assumptions about the model, such as how supersymmetry is 
broken. Virtually all of the existing work on supersymmetric cold dark matter has been done in the context 
of minimal supergravity, mSUGRA, which reduces the hundred unknown parameters to four and one sign. 

But even in the extremely restricted framework of mSUGRA, it may not be possible to determine the relic 
density in general; and there exists very little research on general mSUGRA models. Usually, the mSUGRA 
parameter space is further restricted to those narrow regions in which the LSP is all of the dark matter — i. e. 
the thermally-produced contribution to the relic density of the LSP (denoted \) 1S such that f2 x = f2 c dm- 
However, even in this extremely constrained subset of all possible, extremely simplified (mSUGRA) super- 
symmetric models, there are regions in which the computation of the relic density of the LSP is not known to 
be possible. For example, mSUGRA models in the 'higgs pole' region with large tan j3 suffer from enormous 
computational difficulties |9(. 

Perhaps (not too) surprisingly, we have found the relic density calculated within the framework of 
mSUGRA not to be representative of supersymmetric models even slightly more general than mSUGRA. 
Contrary to assumptions common in the literature, the mSUGRA results are not even indicative of models 
where the LSP is almost entirely Bino 4 . For example, it is argued that in mSUGRA the masses of squarks 
are usually irrelevant to the calculation of the relic density, but that a small change in the mass of the LSP 
can have a dramatic effect fl. This suggests that the relic density of X crucially depends on m v and not so 
much on the other parameters of the model. 

However, when the assumptions of mSUGRA are relaxed even slightly, the correlation between m x and 



4 The Bino is the superpartner of the standard model U(l) gauge boson — a superposition of the photon and Z. A characteristic 
of mSUGRA models is the tendency for an almost purely-Bino LSP. 
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TABLE I: Examples of minimally supersymmetric standard models, consistent with experimental constraints, in 
which the (neutralino) LSP is the wimp dark matter particle candidate. 





fi (GcV) m (GeV) m 1/2 (GeV) m A (GeV) tan/3 


At/mo 


A b /m 


Model l.a 


308.2 413.3 327.3 380.2 11.04 


-0.9846 


-0.9842 


Model l.b 


-441.5 405.5 326.3 337.6 15.07 


0.8234 


0.0224 



the relic density virtually disappears, as many more parameters become important. Figure^ shows the 
relic density versus LSP mass for a few thousand randomly generated, experimentally constrained 5 mini- 
mally supersymmetric standard models which have a Bino-like LSP 6 (no presumptions were made about 
the supersymmetry breaking scenario). The mass of the lightest squark in each model is indicated by color. 
Notice that there is almost no correlation between the LSP mass and relic density. Furthermore, the lack 
of correlation is unchanged even if the mass of the lightest squark were also known. Apparently many more 
parameters must be known to determine il x h 2 . 

Throughout this work, we will make use of a collection of several thousand randomly-generated, con- 
strained MSSMs. Each was generated using the DarkSUSY code [lo|. and are consistent with all known 
experimental constraints 7 . Each model is given in terms of seven (high-scale) parameters, which are: 

• /i, the higgsino mass parameter; 

• m i/2) the (assumed universal, high-scale) gaugino mass parameter; 

• mA, the mass of the CP-odd higgs particle; 

• tan/3, the ratio of the vacuum expectation values of the two higgs doublets; 

• mo, the (assumed universal, high-scale) squark mass parameter; 

• A},, the bottom-sector trilinear coupling; 

• A t , the top-sector trilinear coupling. 

Although these models are extremely simplified MSSMs (having reduced the more than one hundred pa- 
rameters to seven), they do not rely on any particular supersymmetry breaking mechanism and are less 
constrained than mSUGRA models. 

Presumably, calculations done within mSUGRA are not very likely to precisely reflect reality (unless the 
real world were in fact mSUGRA, which is unlikely, and also would be very hard to demonstrate indeed). 
Therefore, it will be necessary to work within a more general framework than that of mSUGRA. But we 
have found that MSSMs even slightly more general than mSUGRA indicate that it may not be possible to 
determine the relic density of the LSP for a general MSSM using data from hadron colliders. 

Consider for example the two (mSUGRA-like) MSSMs listed in Table |I| In Table HU we list the particle 



By 'constrained,' we mean that all of the models are allowed under all current experimental constraints on supersymmetry and 
so are physically realizable. We must emphasize that this may not be the same as frequently-used definitions of 'constrained' 
susy models in the literature. Many authors, for example, consider 'constrained' to refer to exclusively mSUGRA models. In 
addition, some authors further add the restriction that Q x = n c( j m . 
By 'Bino-like,' we mean that the Bino component of the LSP is > 97%. 

These include limits on b — » s-y and (g — 2)u (the early summer of 2004 bounds), together with the requirement that the 
thermal contribution to the relic density be no more than that of all the dark matter. We have not tuned the parameters to 
stringently satisfy all the constraints on (m^, &etc), since these should not affect the general results. 
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TABLE II: Particle spectra for example susy Models l.a & l.b. Notice that although these models differ in relic 
density by a factor of 50, their signals at the LHC are virtually identical. All masses are in GeV. 





n x h 2 


m x 


ntNLSP 


m t - 


(ma - m t -) 


m h 


ran 


m-g gauj 


;ino fraction 


Model l.a 


0.09905 


158 


263 


351 


58.7 


112 


381 


1138 


0.935 


Model l.b 


0.00172 


163 


311 


360 


42.2 


110 


338 


1134 


0.986 



TABLE III: Direct and indirect detection data for example susy Models l.a & l.b. Notice that although these models 
have virtually identical particle spectra, their direct detection signals differ dramatically. 





Px 


Ge signal 6 


Nal signal 


Solar v flux' 2 


Model l.a 


0.2971 


1.18 x 10~ 5 


4.07 x 10" 5 


7.67 x 10" 1 


Model l.b 


0.0052 


3.14 x 10~ 8 


1.07 x 10" 7 


9.13 x 10~ 4 



"We assume the relic and local densities to be proportional. 

'As observed in a 20 keV bin, in units of cpd/kg-keV. An isothermal halo model is assumed. 
c As observed in a 5 keV bin, in units of cpd/kg-keV. An isothermal halo model is assumed, 
''in units of muons/yr-km 2 . 

spectra for each model together with the (thermal) relic density of each. The two models differ in relic den- 
sity by a factor of more than 50: the LSP of Model l.a would account for all of the dark matter while that 
of Model l.b would account for < 2% 8 . Despite the vast difference in their relic densities, the models have 
nearly identical spectra. These represent a proof by example that it will be very difficult if not impossible to 
compute the relic density using optimistic data at hadron colliders even for an MSSM slightly more general 
than mSUGRA. And in practice, this would need to be done for MSSMs much more general than those we 
studied. 

Although the models have strikingly similar spectra at hadron colliders, they give vastly different dark 
matter detection signals. In Table HTT1 we list representative data that would be seen in direct and indirect 
dark matter detection experiments. Notice that the difference in density is quite apparent in their relative 
signal strengths. As described later in this work, the direct detection signals are roughly a measure of the 
product of the local density and the scattering cross section. It is natural that models with similar particle 
spectra would have similar scattering cross sections; such models should be differentiable easily in (and 
perhaps only in) dark matter detection experiments. Based on our results above, we expect that it is not 
possible to calculate the thermally-produced contribution to the relic density of a generic wimp using data 
from hadron colliders — even within the framework of a specific wimp model (e.g. supersymmetry) . Recall, 
however, that there were two general shortcomings of this approach — what we have just described is the 
practical difficulty of obtaining enough of the parameters to make the computation. Even if it were possible, 
however, there remain serious shortcomings of principal. 

2. The Non-Thermal (and Non-Standard) Production of Wimps in the Early Universe 

Thermal production in the early universe is the simplest possible mechanism to explain today's abundance 
of dark matter; it is by no means the only mechanism. There is no reason to suspect that all of the dark 
matter was produced thermally, and any non-thermal contribution to dark matter would yield thermal pro- 



Contrary to the intuition from mSUGRA, notice that the model with the smaller gaugino fraction has the larger relic density. 
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duction estimates to be not only misleading but quite inadequate. The cosmological data which indicates 
today's abundance of dark matter is indifferent to the mechanism by which dark matter particles were pro- 
duced, and there is no known way to determine what fraction of the dark matter was produced thermally. 

Perhaps the most commonly considered form of non-thermal production is the decay of heavy particles 
or condensates into wimps. There are many known examples of non-thermal production mechanisms for 
supersymmetric (LSP) dark matter. For example, inspired by anomaly-mediated supersymmetry breaking 
(AMSB) scenarios, heavy gravitinos produced during reheating could decay into LSPs, creating a cosmolog- 
ically significant LSP relic density [rj. Alternatively, moduli fields (flat directions, e.g. from string theory) 
which acquire mass from supersymmetry breaking could decay into large amounts of LSP dark matter be- 
fore big-bang nucleosynthesis |l2(. Other mechanisms include the late-time decay of Q-balls produced in 
e.g. Affleck-Dine baryogenesis (see e.g. Ref. llJl . or the decay of cosmic strings formed via spontaneous 
breaking of high-scale U(l) gauge symmetries |l4| . These mechanisms could be very important because in 
many viable supersymmetric models (e.g. AMSB), thermal production of the LSP would be insufficient to 
account for the dark matter in the universe alone. 

In fact, several non-thermal production mechanisms, including presently unknown ones, could have oc- 
curred during our cosmological history. And virtually all of these mechanisms critically rely on high-scale 
physics well beyond the reach of foreseeable colliders. Therefore, these mechanisms cannot be excluded 
using collider data and will always exist as possible alternatives to thermal production. Hence, there is no 
reason to suspect that the thermally-produced contribution to the relic density of wimps is the same as the 
actual relic density of dark matter measured cosmologically. These credible objections to thermal freeze-out 
calculations cannot be answered with colliders only sensitive to the effective, low-scale physics. 

But the questions of non-thermal wimp production are by no means the only objections to thermal cal- 
culations. As it turns out, even thermal production can be very sensitive to high-scale physics well beyond 
the reach of colliders. For example, scalar-tensor quintessence dark energy could dramatically alter the relic 
density resulting from thermal production. As shown in Ref. the modified expansion rate of the early 
universe in quintessence models can enhance the relic abundance of wimp dark matter by as much as three 
orders of magnitude. Therefore, high-scale physics could have an enormous effect even on the amount of 
dark matter that is produced thermally in the early universe. 

Therefore, if we hope to account for all the dark matter in the universe we cannot rely on calculations of 
thermal freeze-out. There is no known way to eliminate the possibility of non-thermal production or high- 
scale effects which cannot be accounted for using colliders in the foreseeable future. All of the arguments 
above indicate that it is not possible to definitively determine the relic density of dark matter particles using 
data from hadron colliders alone. At best, with the most optimistic data and less than general models, one 
could hope to compute the 'canonical' thermally-produced component of the cosmological relic density. But 
because of the effects of quintessence, even this calculation cannot be considered conclusive. 

The only way to meet these challenges is to measure the actual local density of wimps and compare this 
with the known local density of dark matter in the halo. This of course cannot be done at colliders. If 
the local density of wimps could be determined, then dark matter could be accounted for in spite of all the 
ambiguities associated with non-thermal production and high-scale physics. Nevertheless, though colliders 
cannot do the job alone, they will be essential to the determination of the actual relic density, as we describe 
below. 

It has not escaped the authors' attention that if it is determined that all of the dark matter in the halo 
is composed of the wimps produced in colliders, then non-thermal production mechanisms could be enor- 
mously constrained by making specific (although perhaps unjustified) correlations between the local and 
relic densities of wimps. We will discuss this further in Section C.l. 



B. Directly Interacting with Wimps in the Local Halo 
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As we have seen, the dark matter problem can be partially addressed with colliders only after overcoming 
enormous obstacles of practicality and ignoring critical objections of principal. At best, colliders alone can 
be used to calculate the cosmological amount of dark matter wimps that may have been produced thermally 
in the early universe. This should be strongly contrasted with dark matter direct detection experiments, 
which hope to measure the actual local density of dark matter wimps in our galactic halo today. 

First proposed by Goodman and Witten [16(, dark matter wimps in the halo could be observed directly 
through their occasional scattering with ordinary matter in earth-based detectors. Extremely sensitive 
calorimeters record the small amount of recoil energy deposited by wimp scattering. Because wimps are 
expected to have masses on the order of a hundred GeV and move relatively slowly in the halo — on the order 
of a few hundred km/s — they would deposit recoil energies up to ~ 200 keV. Direct detection experiments 
measure this recoil, and record the wimp-nucleus scattering rate as a function of recoil energy and time. 

The scattering rate observed in these experiments is a measure of the product of the local wimp density 
and the wimp-nucleon scattering cross section. If one knew the scattering cross section, then these signals 
directly determine the local density of particles. And because we know from astrophysical studies that the 
local density of dark matter is roughly p c dm ~ 0.3 GeV/cm 3 , the local density of wimps directly determines 
the fraction of the local dark matter they compose. It must be noted, however, that there are some important 
caveats involved in this line of reasoning, and these will be discussed below in Section C.l 

But because the scattering cross section is required to compute the local density from direct detection 
data, colliders will be absolutely critical. Only using data from colliders can the cross section be calculated 
to determine the local density. We will illustrate that this could be done through example and explicit 
calculations. 



1. Why the Local Density of Wimps Cannot be Determined with In/Direct Detection Alone 

Let us imagine that a weakly interacting massive particle \ nas been unambiguously observed in dark 
matter detection experiments. Because the interaction rate of wimps with a detector is a measure of the 
local density multiplied by the scattering cross section, low-density wimps with a high cross section are 
indistinguishable from high-density wimps with a low cross section. But it is precisely these two extremes 
that are produced by 'thermal frcczc-out' in the early universe: if the cross section is large then more wimps 
would have annihilated in the early universe and resulting density would be small; alternatively, if the cross 
section is small then freeze-out would have occurred very early and the density today would be high. This 
crude argument suggests that even a very small component of dark matter may have a detectable signal: it 
would likely have a relatively large cross section. This well-known result has been referred to as the 'no-lose 
theorem' in recent work and conferences. Indeed, experimentalists may not lose out on discovering even a 
tiny fraction of the dark matter But 'no-lose'^win. 

To see how little the density of wimps is correlated to our ability to detect them, we explored randomly- 
generated (physically allowable) supersymmetric dark matter models. In Figure [3 we plot the local density 9 
against a representative direct detection signal for some six thousand constrained minimally supersymmetric 
standard models (without assuming any supersymmetry breaking scenario). The lack of correlation is imme- 



This data was generated assuming that the local density is proportional to the cosmological density. 
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Direct Detection Signals 
versus Local Density of Wimps 




£ I , — , . — — . ™j , — „ , l 

£ 10" lxlO" 5 lxl 0" 4 10"' 

Direct Detection Signal in Ge at 20 keV (cpd/kg-keV) 

FIG. 2: The local halo densities of several thousand constrained MSSMs as a function of direct detection signal 
strength (in germanium). To generate this plot, the galactic and cosmological fractions of wimp dark matter were 
assumed to be same — see Section C.l for a discussion of the implicit caveats of this line of reasoning. Experiments 
currently in planning or under construction may be able to observe signals of the order 10 -4 cpd/kg-keV. Notice that 
any particular signal strength is consistent with a large range of local densities. 

diately obvious. Notice that for any particular signal strength the local density fluctuates over at least two 
orders of magnitude. Hence, in accordance with the no-lose theorem, the discovery of wimps in the galactic 
halo tells us very little about how much of the local dark matter they compose. Figure [21 indicates that a 
wimp discovery in the near future could easily represent a negligible fraction of the local dark matter, even 
1% or less. 

Because direct detection signals are universally proportional to the local wimp density, it is clear that no 
combination of direct detection data could resolve the density — there would be no data independent of this 
parameter. What is not as obvious, however, is that indirect detection data may not resolve the density 
cither. There are two principal methods of indirectly detecting wimps in the halo: the study of high energy 
neutrinos from the annihilation of wimps in the sun or earth and the study of cosmic rays. The flux of 
high energy neutrinos from wimp annihilation in the sun or earth offers little more insight into the local 
density of wimps because — considering the age of the solar system — the wimp annihilation rate is controlled 
by the wimp-capture rate (is| : this is essentially a measure of the elastic scattering rate, and so cannot be 
used to separate the local wimp density from the scattering cross section any more than direct detection 
experiments. 

The study of the galactic gamma-ray spectrum or the anti-matter excess in cosmic rays — although not 
directly related to scattering rates — cannot be used to determine the density of a wimp for different reasons. 
Specifically, the total flux of particles from wimp annihilation in the halo is proportional to a local, unknown 
'boost factor' 10 which can enhance the flux of antiparticles and gamma rays by one or two orders of magni- 
tude (for a recent review see Ref. [l9|). Therefore, these indirect searches are sensitive to the shape of the 
wimp density distribution, but insensitive to the actual local density of wimps. 

From our discussion above, it should not be surprising if the local density of wimps cannot be determined 



The 'boost factor' is related to galactic confinement times and other astrophysical issues. It collectively represents many 
(typically local) astrophysical uncertainties. 
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TABLE IV: Examples of minimally supersymmetric standard models, consistent with all current constraints, in which 
the LSP is the wimp dark matter particle candidate. 





fj, (GeV) m (GeV) m 1/2 (GcV) m A (GeV) tan/3 


At/mo 


A b /m 


Model 2.a 


463.4 593.8 371.1 428.2 18.39 


-0.1844 


0.4366 


Model 2.b 


-412.1 435.5 372.7 337.4 23.33 


-0.0119 


0.0707 



TABLE V: Direct and indirect detection data for example susy Models 2. a & 2.b. Notice that the two models have 
nearly identical direct and indirect detection signals although their relic densities differ significantly. 





„ o 

Px 


m x (GeV) 


Ge signal 6 


Nal signal 


Solar v nux d 


Galactic p flux e 


m fi (GeV) 


Model 2. a 


0.2856 


184 


3.50 x 10~ 6 


1.18 x 10" 5 


0.121 


4.58 x 10~ 10 


591.3 


Model 2.b 


0.1044 


185 


4.09 x 10~ 6 


1.38 x 10~ 5 


0.106 


9.92 x 10~ 10 


432.2 



"We assume the relic and local densities to be proportional. 

'As observed in a 20 keV bin, in units of cpd/kg-keV. An isothermal halo model is assumed. 
c As observed in a 5 keV bin, in units of cpd/kg-keV. An isothermal halo model is assumed, 
''in units of muons/yr-km 2 . 

e Measured at 3 GeV, in units of p/cm 2 -sr-s-GeV. 

in general from any collection of dark matter detection experiments. Indeed, we have found many examples 
where models with quite different relic densities show nearly identical signals for every type of direct and 
indirect detection experiments. Consider the two (mSUGRA-like) supersymmetric models listed in Table HVl 
We calculated the in/direct detection signals and relic densities for each model using the DarkSUSY code 
[lfj : these are listed in Table [V] Notice that although the two models give rise to virtually indistinguishable 
in/direct detection data, they differ substantially in relic density. 

Also listed in Table [V] is the mass of the scalar up-squark for each model. Although it could be mea- 
sured only at colliders, we list rriu among these experiments for a purpose that will be better understood 
later in this work. As we will see, the masses of the squarks are critically related to the scattering cross 
section. Heuristically, scattering is typically dominated by squark exchange and is roughly proportional to 
l/(m| — m^.). Therefore, two models with the same detection signals will have densities such that 

(£xh K~)i - K)i m 
( Px h ~ K) 2 -K) 2 - 1 ) 

This back-of-the-envelope analysis works remarkably well for the two models, as can be easily verified nu- 
merically. 

Ours and similar studies suggest that no combination of dark matter in/direct detection data can be 
used to determine the local density of wimps alone. Therefore, it is critical that wimps also be observed in 
colliders. By producing wimps at colliders and studying their particle physics, the wimp-nucleon scattering 
cross sections may be calculated. Combining the scattering cross sections with direct detection data would 
immediately give the local density. 



C. Solving the (Local) Dark Matter Problem 
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Our studies suggest that the dark matter problem can only be solved in general if wimps are observed 
and studied at both colliders and direct detection experiments. In Section A, we showed that colliders could 
at best estimate the (quintessence-free) thermally-produced component of the cosmological relic density of 
wimps. But these analyses require enormous lack of generality through oversimplification (such as assuming 
mSUGRA in supersymmetric models) and rely on high-scale physics which could not be observed at colliders 
(such as the non-existence of quintessence or additional high-scale gauge symmetries). In contrast, direct 
detection experiments are sensitive to the actual local density of dark matter wimps. But these experiments 
cannot measure the local wimp density without first knowing the relevant scattering cross sections — and this 
can only be determined using collider data. 

To determine the cosmological significance of wimps observed in both direct detection experiments and 
at colliders, we must measure their density in the galactic halo. This can be done by combining scattering 
cross sections determined using collider data with data direct detection experiments. But first it must be 
established that the particles observed in the two experiments are in fact the same 11 . There is only one 
known way to suggest that the two signals correspond to the same wimp: by showing that each has the same 
mass. Fortunately, it may be possible to determine the mass of particles observed in each type of experiment 
independently. 

Perhaps surprisingly, it appears much easier to determine the mass of wimps using direct detection ex- 
periments than with colliders. This is because there exist model-independent ways of determining the wimp 
mass with direct detection data alone; these methods are described in Section III.B. At least one of these 
methods, such as that presented herein, could always be used to determine the mass (in principal). In strong 
contrast, there exists no general procedure to determine the mass of wimps produced at hadron colliders. 
This is because, being weakly interacting and stable, wimps indicate their presence as missing energy only: 
mass differences are easily measurable, but masses themselves require intensive statistical analysis. Such 
statistical techniques have been developed to determine the mass of the LSP for some types of mSUGRA 
models (see e.g. Ref. 0|). However, there is no known way to measure the LSP mass in the general MSSM. 
It may not even be possible to determine the mass of an arbitrary wimp produced at colliders in a model- 
independent way. Even in the (albeit extremely favorable) situations of some supersymmetric models where 
kinematics alone allows one to determine the LSP mass, the methods require enormous amounts of data. 
Therefore, it could be many years before there is a determination of a wimp mass from hadron colliders 12 . 
But such a determination is crucial to the identification of the wimp seen in direct detection experiments. 

After the particle has been identified and its mass is known, then collider data can be used to estimate 
(or bound) the wimp-nucleon scattering cross section; this allows one to determine the actual local density 
of wimps. In contrast to the 'thermal' relic abundance calculations described above, it is found that only 
minimal data from colliders is required to absolutely bound the scattering cross section for general (super- 
symmetric) models. Therefore, even sparse collider data could immediately translate into strong bounds on 
the actual local wimp density. 

We will illustrate how these estimates are made in the framework of a supersymmetric (neutralino) wimp 



11 Furthermore, it must be established that direct detection experiments observe only a single type of wimp particle: maybe 
multiple wimp components exist in the halo. We are unaware of any concrete way of establishing this fact, other than perhaps 
observing consistency in the mass calculations described in this work. 

12 Of course, soon after signals of new physics are reported at hadron colliders, such mass determinations will be reported; 
but they will likely depend on many restrictive, unconfirmed assumptions (such as mSUGRA). But solving the dark matter 
problem is too important to risk these types of misguided errors. 
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for the most general softly-broken MSSM. We will see that not only can the scattering cross section be 
bounded strongly using even minimal data, but also that this bound rapidly approaches the true value as 
more collider data is obtained. These conclusions are true without relying on any assumptions about super- 
symmetry breaking, relative sparticle masses, or other restrictions on the MSSM. 

Combining estimates or bounds of the wimp-nucleon scattering cross sections from colliders with direct 
detection data immediately allows one to estimate or bound the actual local density of wimps. Up to some 
caveats described below, this 'measured' local density of wimps would allow one to determine what fraction 
of cold dark matter they represent, thereby addressing their cosmological significance. To the extent that it 
is determined that p x , measured with colliders and direct detection experiments, is equal to p c d m , measured 
by the orbital velocity of nearby stars in the Milky Way, the dark matter problem will be solved. 

1. Caveats and Ambiguities About the Local Density of Wimps 

Because direct detection experiments are sensitive to local, small-scale structure in the galactic halo, our 
existing knowledge of the ambient halo density may not be sufficient to solve the dark matter problem in 
general. Our knowledge of the local halo density — that is, the fact that /d c dm ~ 0.3 GeV/cm 3 — is based on 
relatively large-scale surveys of star velocities in the Milky Way; these measurements are not very sensitive 
to small-scale structure in the halo such as local density perturbations. 

Most of the small-scale structure considered in the literature involves local, high-density regions of dark 
matter within the halo. Although these structures would make it easier to discover dark matter, they make 
it nearly impossible to assess what fraction of the halo is composed of a discovered wimp. (But of course, 
we could be in a low-density region.) 

There are several types of small-scale halo structure which may effect direct detection experiments. For 
example, some authors have proposed that the halo may be clumpy or contain caustic structures such as 
rings or shells (see e.g. Refs. |2l|,|22j]). Alternatively, the earth may be within a stream of dark matter. This 
situation has been suggested by studies of the Sgr A galactic tidal stream; it has been estimated that this 
stream may increase the local halo density by 0.3 — 23% relative to the ambient density 23]. These small- 
scale perturbations in the dark matter density could have significant effects on direct detection rates — and 
nearly fatal effects on our ability to resolve the composition of dark matter. 

Fortunately, there may exist ways to eventually exclude or compensate for most of these perturbations. 
For example, dark matter streams or caustics would produce anisotropy in the local wimp velocity profile; 
these may be identified using directional detection experiments like DRIFT Q . A clumpy dark matter halo 
may possibly be identified by studying the time-dependence of a wimp signal 13 . These ambiguities will need 
to be addressed before the dark matter problem has been put to rest. Fortunately, because these obstacles 
could possibly be surmounted in the future, they are much better than those associated with non-thermal 
(and non-standard) production in the early universe. A determination of the actual local density of dark 
matter particles is as close to a cosmological result as we can get. 

Unfortunately, there does not exist any well-understood or accepted physical argument which allows one 
to go from the the local halo density to the cosmological relic density or vice versa. The relationship between 
the two involves many details of galaxy formation and structure that remain to be understood. Recall 



It may not be possible to observe fluctuations in the scattering rate caused by a clumpy halo, considering the time scales 
which may be involved. These types of questions will not have definite answers until more work has been done to understand 
the physics of realistic clumpy halo models. 
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that the two principle sources of evidence for the existence of dark matter, cosmology and astrophysics, 
are both independent and distinct: it is a hopeful assumption that the dark matter required by the cosmic 
microwave background and large scale structure formation is precisely (and entirely) the same material that 
is required by the dynamics of galaxies. Therefore, even if all the ambiguities about the local halo structure 
were resolved and all the local dark matter was identified, it not obvious that all of the dark matter in the 
universe has been accounted for. However, these ambiguities would not prevent our complete understanding 
of the composition of the Milky Way's dark matter halo. Of course, we do expect that cosmology and 
astrophysics are measuring the same dark matter, but the issue is so important that this should not be 
automatically assumed. 

2. Constraining the High-Scale Physics of Non-Thermal Wimp Production 

Recall that there exists many non-thermal mechanisms to produce wimps in the early universe or enhance 
the number of wimps thermally produced. If and when the local density of wimps is established by the 
procedure described above, we would be tempted to infer the relic density by associating the relic density of 
dark matter, f2 c d m /i 2 ~ 0.1, with the local density of dark matter, p c( jm ~ 0.3 GeV/cm 3 . //this assumption 
were made, then up to the caveats described in the previous section, one could infer the actual relic density 
of wimps. As described above, this is not a trivial assumption. 

And if there existed sufficient data to approximate the amount of wimps produced thermally in the early 
universe, this could be compared to their inferred relic density. Such analyses could enormously constrain 
the existence of many forms of high-scale physics. For example, an agreement between the two results 
would imply that all of the dark matter was indeed produced thermally and without enhancement from 
quintessence. Alternatively, any discrepancy between the two would imply the need for new physics at the 
high-scale. This type of analysis would be a unique window to high-scale physics. 

Therefore, after wimps are discovered and their local density determined, thermal freeze-out calculations 
may improve our understanding of the very same physics which today prevents freeze-out calculations from 
helping us understand dark matter. 

III. WHEN WIMPS ARE FIRST DISCOVERED: MAKING THE MOST OF DIRECT 

DETECTION DATA 

As we have shown, direct detection data is critical to solving the dark matter problem. When data first 
becomes available, it will be necessary to determine the wimp mass and separate the various independent 
contributions to the scattering rate. To elucidate this procedure, we will briefly review the standard formalism 
for the direct detection signal. This will allow us to demonstrate methods to optimize data, and to introduce 
an explicit technique to determine the wimp mass. This discussion will also introduce the formalism used 
later in this work to bound the LSP-nucleon scattering cross section. 



A. The Scattering Rate of Halo Wimps with a Solid-State Detector 
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Let us imagine that the wimp \ has been observed in direct detection experiments. It will be helpful for 
us to review the explicit form of the differential interaction rate for a particular detector measured at recoil 
energy q. Let the detector in question be composed of nuclei labeled by the index j, each with mass fraction 
Cj. Then, the differential rate of wimp scattering at recoil energy q is 14 , 



dR 
dQ 



Q = q Jv^iq) 



^-dv |f? (?) [Zjfp+iAj - Z i)fn\ 2 + pj^+i) [ a i^«o (l)+ a o S 3ii {q)+aia S Jol (q)] 

(2) 



where v m i nj (q) is the minimum velocity kinematically capable of depositing energy q into the j th nucleus; 
f(v,t) is the local velocity distribution function for wimps in the galactic halo; Fj(q) and Sj mn (q) are 
nuclear form factor functions for coherent and incoherent scattering, respectively; Zj and Aj are atomic and 
mass numbers; Jj is the nuclear spin; a\ = a p + a n and ao = a p — a n ] and the constant parameters f p>n 
and a Pt n describe the coherent and incoherent wimp-nucleon scattering cross sections, respectively 15 . The 
interaction parameters can be calculated for any explicit wimp model, regardless of the type of wimp under 
consideration. 

In general, this expression for the scattering rate is a second order polynomial in the four (unknown) 
interaction parameters / Pi „ and a PiTl . To highlight this, it can be recast in the suggestive form, 
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(3) 

It is clear from the expression above that using data from 

1. different detector materials (varying the mass fractions, nuclear form factors, nuclear spins, and mini- 
mum velocities) and 

2. different recoil energies (varying the nuclear form factors and minimum velocities), 

one can invert equation (|3J) to solve for ^fp^f p ^ n and yfp^ P .n if the halo velocity distribution and m x were 
known. That is, given a halo model and wimp mass, one can use data from different detector materials and 
different recoil energies to determine y/p^f p ,n and ^fp^a p ^ n (up to quadratic ambiguities). 

Although both the coherent and incoherent (scaled) interaction parameters are in principle calculable, it 
may not be possible to determine them all when wimps are first discovered. It is perhaps even likely that 
there will not exist enough data initially to solve for both. For example, if wimps are observed with detectors 
largely insensitive to spin effects — as is the case for 76 Ge-enriched detectors — then it will not be possible to 
determine the spin-dependent contribution to scattering. 



14 For a detailed account of equation J2J see e.g. Ref. 

15 Imprecisely, coherent scattering is sometimes called 'spin-independent' and incoherent scattering 'spin-dependent.' 
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The Difference Between Summer and Winter 
Differential Wimp Interaction Rates 
For Various Energy Bins 

3x10-" T — 




Energy Bin (keV) 

FIG. 3: The difference between direct detection signals in June and December as a function of recoil energy for several 
detector types. This plot was generated for a model with m x ~ 161 GeV. 

Notice that in equation Jjyl, the coefficients of a p and a n are linearly independent combinations of the 
incoherent form factor functions Sj kl (q) ; this is in contrast to f p and f n which are each multiplied by the 
same form factor function Fj(q). This implies that rates observed at different energies can be used to 
separate a p from a n with any single detector capable of measuring incoherent scattering, but rates observed 
at different energies cannot in general distinguish f p from /„ . It is apparent that to isolate coherent scattering 
contributions from neutrons and protons one must use data from detectors made of distinct materials. 

It has not escaped the attention of the authors that knowing the interaction parameters up to an unknown 
factor — yJ~Px — can Si ye enormous insight into the particle physics of \. For example, if \ is the LSP then 
the ratios a p /a n and a p /f n could lead to important constraints on tan/5, the degeneracy of squark masses 
and mixing parameters, and perhaps contain other information as well. This could be very important for 
interpreting collider physics and disentangling the MSSM. 

B. Determining the Mass of the Wimp 

For a long time there was only one known way to determine the mass of a wimp using direct detection 
experiments. As the earth orbits the sun, its velocity through the galactic dark matter halo varies between 
roughly 190 and 250 km/s |8|. This causes annual modulation in the scattering rate. It has long been known 
that the amplitude of this annual modulation can change sign between high and low recoil energies; if so, 
there is a particular energy, the 'crossing energy,' at which no annual modulation is observed. This can 
be seen in Figure for several detector materials. In a review article in 1988, Ref. [2^. Primack et. al. 
showed that the crossing energy is directly related to the wimp mass 16 . The method is generally robust if 
the crossing energy is in fact observed. 



Although it seems unlikely to have originated in the review article , we ha ve been unable to find any previous author mentioning 
the effect. For recent and more thorough descriptions sec Refs. I2 r ll2(l . 



18 




FIG. 4: The annual modulation crossing energy as a function of wimp mass for several thousand randomly-generated 
constrained MSSMs. Notice that the effect would be difficult to observe for light wimps or with a detector composed 
of light nuclei, e.g. sapphire (AI2O3). Notice the ambiguity which arises in the case of detectors which are composed 
of multiple nuclei. 



Unfortunately, depending on the masses of the wimp and detector nuclei, the crossing energy may not 
be observable in realistic experiments. Also, the effect can be complicated by the presence of multiple nuclei 
types in the detector. If the wimp is light then the crossing energy could be lower than any planned detector's 
threshold: effectively, the annual modulation will peak in June for all of the detector's energy bins (including 
the lowest) 17 . This is also the case for detectors made of relatively light nuclei, e.g. sapphire (AI2O3). In 
Figure 01 we plot the annual modulation crossing energy as a function of mass for detectors made of sodium 
iodide, germanium, and sapphire. Notice that each of the limitations just described is clearly illustrated in 
the figure. 

Recently, we found an alternative technique to determine the mass of a discovered wimp using data from 
direct detection experiments alone 0, Q . This method is based almost entirely on kinematics and will work 
in principle regardless of what particle the wimp is. 

The technique presented below may not be entirely unlike those used implicitly in existing analyses (see 
e.g. Ref. 28|). However, we are unaware of any formal or explicit presentation of these techniques in the 
literature. 



1. Using Kinematical Consistency to Determine the Wimp Mass from Direct Detection Data 

We have found that the wimp mass can be determined if signals are observed in either a single detector 
that is capable and sensitive enough to observe the incoherent contribution to the scattering rate in several 



Notice that this is precisely what has been observed in the DAMA experiment 1271 . Hence, if the DAMA signal is in fact 
a signal, it favors a (surprisingly) low wimp mass. However, only an upper bound is obtained using this observation alone 
because the experiment did not observe an energy bin without annual modulation. 
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Algorithm to Determine the Mass of the Wimp 
Given Direct Detection Rates 
for both Nal and Ge Detectors 

(model is ATLAS_SUSY02) 
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Algorithm's Scan 
True Wimp Mass 
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Mass of the Wimp used for Calculation 



FIG. 5: The function C( m x) where the wimp is the neutralino in the MSSM example specified by ATLAS susy point 
2 [gjj. The method is equally applicable for arbitrary models of wimp dark matter; in particular, LSP dark matter 
models with arbitrary soft parameters. 



energy bins, or in at least two detectors made of distinct materials in several energy bins. This is done as 
follows. As we have shown, combining data from multiple detectors and/or from multiple recoil energy bins 
allows one to isolate the relative scattering contributions from coherent and incoherent interactions with 
neutrons and protons. To do this, however, the mass of the wimp and the velocity profile were required. 
If the halo velocity distribution were known, then only the wimp mass would be required. Let us assume 
this to be the case. Then, given the wimp mass, equation @ can be inverted to solve for ■s/Pxfp,n and/or 
^J~p^a p ^ n . Because many direct detection experiments observe scattering rates for many recoil energies, we 
can expect the system of equations to be over-constrained soon after wimps are first observed. 

But what if m x is not known? The system could still be inverted with an assumed mass m' x ; but if 
m' 7^ m x then distinct sets of data would generally give inconsistent results. That is, because the interaction 
parameters are constants all minimal, linearly independent combinations of measurements used to solve for 
the scaled interaction parameters will agree if the correct mass were used in the derivation. However, if 
an arbitrary m' were used to solve for these parameters, different calculations will not in general agree. 
Therefore, the requirement of consistency in these calculations can be used to identify the correct wimp 
mass 18 . 

This motivates us to define a 'kinematical consistency' function Q(m' x ) which compares the values of 
the scaled interaction parameters obtained using different independent subsets of data as a function of m' x . 



This technique is quite similar to x 2 - m inimization or maximum likelihood analys es. As stated earlier, similar statistical 
techniques are already used by experimental groups to determine m x (see e.g. Ref. l2gl ): but we are not aware of any formal 
presentation of these in the literature. 
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Specifically, if there exists sufficient data sufficient data to solve for ^/p^/p^n 19 , then let ((m' x ) be given by 

C(m' x ) = VPx {(/p(0 - fpti)) 2 + similar terms j , (4) 

where the indices i,j represent distinct minimal sets of data used to compute the constants using the value 
m' x for the wimp mass. It is necessary that C( m x) = when m' x = m x . That is, independent calculations 
must agree if the correct mass is used. We should mention that this may not be a sufficient condition; but 
although we have not found any way to demonstrate that m x is the unique root of £(m^), we have found no 
example where it has multiple roots. 

To determine the wimp mass, one varies m' x until C( TO x) = 0- To test this technique, we applied it to our 
collection of several thousand randomly-generated, constrained MSSMs. The correct mass was determined by 
this method for every single model tested. Figure [S] illustrates a typical plot of £(m^), ignoring experimental 
uncertainties and resolutions. The minimum of C( m ' x ) a ^ the correct mass is far lower than £ evaluated at 
(m x ±5GeV). 

There still remain important questions about how this technique can be used in practice. Specifically, our 
studies have not yet considered the effects of experimental uncertainties, backgrounds, or resolutions. Also, 
we do not yet have a complete understanding of the requirement that the halo model be known. We have 
generally found that for direct detection data simulated using a 'smooth' halo profile, the correct mass is 
obtained even if a different (smooth) profile is used in the analysis. However, the effects of using a smooth 
profile in the analysis if the halo is in fact clumpy or anisotropic is not fully understood. 

IV. THE CASE OF SUPERSYMMETRIC (NEUTRALINO) DARK MATTER 

The most popular and perhaps best-motivated wimp candidate for cold dark matter is the lightest su- 
persymmetric particle (LSP) predicted in supersymmetric extensions of the standard model which conserve 
i?-parity 20 . Indeed, supersymmetric dark matter was predicted well before it was known that non-baryonic 
dark matter was needed to solve the dark matter problem. In many MSSMs allowed by experimental con- 
straints the LSP is the neutralino, which we denote x; it is the supersymmetric partner of the neutral gauge 
and higgs bosons. 

Based on the arguments presented above, if the LSP is observed in direct detection experiments and 
produced at colliders then its actual local density can be determined if its scattering cross section (i.e. its 
interaction parameters) were known. Similarly, the local density of x can be bounded if bounds on the inter- 
action parameters were known. As it turns out, one can bound the parameters using even only sparse data 
from colliders, without making any assumptions about the MSSM lagrangian parameters. In this section, 
we will derive this bound explicitly for the most general MSSM. 



The analysis is of course improved if both the coherent and incoherent (scaled) parameters are determinable. If so, it should 

be obvious how to include this information to strengthen the result. 

For a comprehensive review of supersymmetric dark matter see Ref. llSl . 



A. Neutralino Interaction Parameters 
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Dark matter experiments can be used to determine the effective %-nucleon interaction parameters ^fp^a p ^ n 
and/or J~f>xfp,n- These are of course functions of the %-quark (and x-gluon) scattering cross sections. To 
tree level, the %-quark scattering cross sections are given in terms of the following diagrams. 

X ft ft XX ft X ft 




X ft X ft X ft X ft 

Given a completely specified MSSM lagrangian it is not difficult to evaluate all of these diagrams and 
determine each of the interaction parameters. But these will depend on many of the model parameters; in 
particular, they are functions of 

1. the gauge content of x (its photino, zino, and higgsino components), 

2. many squark masses and mixing angles, 

3. tan P, the ratio of the vacuum expectation values of the two higgs doublets, 

4. higgs mass parameters (only for the coherent scattering). 

It must be emphasized that most of these parameters will be very difficult to measure in practice (especially 
at hadron colliders). Today, there do not exist general, model- independent methods of determining most of 
the required parameters at hadron colliders, which is all that will be available for well over a decade. 

Therefore even years of excellent collider data may be insufficient to properly calculate these scattering 
cross sections without making unjustified assumptions about supersymmetry breaking, for example. Al- 
though it may be tempting to interpret any initial collider data in the limited framework of particular 
supersymmetry breaking scenario, e.g. mSUGRA, this can be very misleading. Such analyses trade param- 
eters to purchase an inherent lack of authority. And the possibility of understanding dark matter is too 
important to oversimplify. 

However, as we will show, it is possible to obtain authoritative limits on the LSP-nucleon scattering cross 
sections without resorting to any model for supersymmetry breaking. More astonishingly, this can be done 
using only extremely limited data from colliders. 

B. An Absolute Bound on the Incoherent Scattering Parameters of the LSP for an Arbitrary 

MSSM (Using Minimal Input from Colliders) 

Using only a lower limit on the mass of the lightest squark and bounds on tan f3, we find a strict, model- 
independent upper bound on the incoherent x-quark interaction parameters a Ut d, s - This translates into a 
strict upper bound on the parameters a Pjn ; together with direct detection data to solve for ^J~p^a p ^ n , this 
bound determines a strict, model-independent lower bound on the local density of neutralino dark matter. 

Recall that functions (and in fact linear combinations) of the x-qu^k interaction parameters 

o<u,d,s 21 - Let us derive a strict upper bound for a u . Although all three x-quark parameters are required to 



The explicit relationships between a P;Il and a u $ s are given in Ref. Il8l . 
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compute a Pt „, it will be obvious how arguments similar to those for a u can be made for a dyS - 

In terms of the lagrangian parameters of an arbitrary softly-broken MSSM, the axial-vector couplings of 
X to the u, d, s-quarks are given by 22 



16m^ Hl H2 8 t-r 1 to? -(m x +m„) 2 ' 

3 2 _ fAr2 ^^1^ (*«ai) a + W,y) a 
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o d j v2mv^sm/i 2 



W dl3 = -gV2\t & nO w Nl{U R ® d ) - . . 
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In the expressions above: TIl.r are 3x6 projection matrices given in the basis (ul, cl, th, ur, cr, £r); 9„ is 
the unitary matrix which diagonalizes Af 2 so that M 2 dlag = O^M^G-u 23 ; the subscript j on qj corresponds 
to the flavor and handedness of the squarks, so that j — I, . . . , 6 corresponds to {ul, cl, tLi ur, cr, tn); and 
gauge content of x is given by X = N S \B) + N W \W) + N &1 \H X ) + N 6a \H 2 ) 2i . 
Expanding the expression for a u , we obtain 
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(5) 



Similar expressions describe the scattering with d, s-quarks. 

Consider the contribution to a u from squark exchange — the terms proportional to 1/ (to?. — (t77 x + m u ) 2 ). 
Generally, this will be dominated by scattering via the exchange of u, d, s— squarks 25 . Suppose that there 
exists a lower bound on the mass of the lightest squark, rriq e — which is of course also a lower bound on 
the mass of each squark. This does not require that the lightest squark has been found; rather, it requires 



These expressions are given in terms of the conventions of Ref. Ilgl : please refer to its appendix for a more detailed account. 
We have chosen the basis for the squark-mass matrices so that M u ^ are diagonal. 

As explained briefly before, the fields B and W are the supersymmetric partners of the standard model U(l) and U(2) neutral 

gauge bosons. They are linear combinations of the supersymmetric partners of the photon and Z. 

This statement is of course sensitive to the distribution of squark masses and mixing parameters (unknown today). 
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that there exist no squark which has a mass of less than rriq e — so existing, negative collider results would be 
sufficient. If this is so, then 

1 1 



< 



m l - ( m x + m lt ) 2 m| - (m x + m u ) 



If the squark masses are nearly degenerate, as is often the case, then the inequality becomes a good approx- 
imation. It must be emphasized that our results do not rely on any assumptions about the relative masses 
of the squarks. However, the limit will be more stringent for models where all the squarks have similar 
masses 26 . Substituting l/(m~ ( — (m x + m u ) 2 ) for l/(m|. — (m x + m„) 2 ) in the summand of equation we 
can greatly simplify and limit the expression. Notice that the only j-dependence remaining in the expression 
comes from the squark mass mixing parameters, {J^L,R^u,d)ij- But here, we recall that unitarity implies 
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Therefore, we obtain the strict upper bound 
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where a= {ot w ) is phase between Ng and iV^ (TV^) 27 . 

This expression has seven real unknowns including sin (3. By the normalization of the neutralino wave 
function, \N § \ 2 + \N W \ 2 + \N 6i \ 2 + \N S J 2 = 1, the six unknowns { | N s | , | \ , | N= ± \ , \ \ , a w , } form 
a compact parameter space. Therefore, a u can be absolutely maximized with respect to all these numbers 
(unknown today). The only unknown which cannot be marginalized is sin/3. To maximize a u with respect 



If the squark masses are widely separated, then the derived upper bound will significantly overestimate the scattering 
parameters. However, the bound could be significantly improved if, for example, the non-degeneracy of squarks was a 
fact known experimentally. 

Note that we have allowed the soft lagrangian parameters to be complex. 
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Lower Bound (and Estimate?) 
of the Local Neutralino Density Using Nal Data 
with Bounds on Lightest Squark Mass and tan/3 




True Local Density (GeV/cm 3 ) 



FIG. 6: This plot compares the lower bound (and estimate?) of the neutralino local density computed using the strong 
upper bound for a p , n to the true density for thousands of models. Notice that the procedure correctly determined a 
lower bound for the local density for every model. Unfortunately, sometimes the lower bound was well below the true 
density; this is the case when the lower bound on the lightest squark mass was poor. Recall that the local density of 
halo dark matter is p c d m ~ 0.3 GeV/cm 3 . 

to sin (3, a lower bound of sin j3 must exist — corresponding to a lower bound of tan /3 28 . Using this bound in 
place of sin/? in the expression above, we obtain an absolute upper bound for a u . 

Notice that in our work above we allowed all of the gauge content of the neutralino to be unknown and 
admitted possibly complex, arbitrary soft parameters and phases. In short, the analysis applies to the most 
general softly-broken super symmetric standard model. 

It is important to note the flexibility and improvability of the derivations for these bounds. If, for example, 
the masses of several light squarks were known then one could greatly improve the bounds by including these 
in the explicit expression for a q and maximizing it relative to the remaining unknown parameters. In this 
manner almost any additional knowledge can be added to arrive at stronger results. Therefore, not only do 
these bounds grow more restrictive with increasing knowledge of the MSSM, but they continue to approach 
a realistic estimate of the interaction parameters. It should also be noted that the approach is easily 
generalizable to the case of a (neutralino) metastable next-to-lightest supersymmetric particle (NLSP) if it 
is the dark matter observed in these experiments. 

The coherent parameters are similar to the incoherent ones except that they also contain higgs exchange 
at tree-level. This implies that in addition to squark masses, mixing angles, tan/?, and gauge content of 
X, one must know the higgs masses. Therefore, in general, less is required to compute a Pj „ than / Pjn . But 
even with allowing the higgs masses to be 'known,' we have not been able to find robust bounds for the 
parameters / Pi „ comparable to those for a p ^ n . It appears that much more data from colliders will be required 
to compute or bound / p , n . 



Both the upper and lower bounds of tan /3 are important. This is because a lower bound of cos /3 — effectively an upper bound 
on tan j3 — is needed for a<j s . 
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C. A Strong Lower Bound on the Local Density of Neutralino Dark Matter 

In the previous section we showed that even weak bounds on the mass of the lightest squark and the 
range of tan/3 allows one to obtain strict, model-independent upper bounds on the incoherent scattering 
parameters of x f° r neutrons and protons. Because direct detection data measures the product of the 
scattering parameters and the local density of ncutralinos in the halo, the upper bound on a p , n immediately 
translates into a strong, model-independent lower bound on the local neutralino density. 

To test this, we considered several thousand randomly-generated, experimentally constrained MSSMs. For 
each of these models, we calculated the interaction rates for a Nal detector in twelve recoil energy bins. This 
(idealized) data was used to compute the mass of the LSP, using the kinematical consistency function, and 
to solve for yfp^a p , n . Upper bounds were calculated for a p . n assuming 10% uncertainty in tan/3 and a lower 
bound on the lowest squark mass of either 200 GeV or the actual mass of the lightest squark, whichever 
was less. The specific gauge content of the neutralino was taken to be known for each model for the sake 
of computational simplicity 29 . Using the upper bounds for a Pi „, we obtained a lower bound on the local 
density of neutralino dark matter, p x . 

Figure illustrates the results of using this algorithm for each of the randomly generated MSSMs. Notice 
that the estimated local density is always strictly less than the true local density — as required by it being a 
lower bound. Also, notice that for many models the lower bound was not such a poor estimate. This will be 
the case, for example, when the lightest squark mass is near or below the (assumed to be 200 GeV) lower 
bound and the squarks are relatively degenerate. 

V. CONCLUSIONS AND DISCUSSION 

We have seen from general arguments and robust specific examples that a discovery of dark matter 
particles in the galactic halo or at colliders cannot by itself solve the dark matter problem. To determine 
the cosmological significance of discovered wimps, their local or relic density must be determined. We have 
proved by example that this cannot be done using data from colliders alone or in/direct detection experiments 
alone. However, we have shown that it may be possible to determine the actual local density of discovered 
wimps by combining data from both colliders and detectors: detectors to measure the product of the local 
density and the scattering cross section, and colliders to estimate or bound the scattering cross section. Up 
to specific and possibly resolvable astrophysical uncertainties regarding the homogeneity and isotropy of the 
local dark matter halo, this may be sufficient to solve the local dark matter problem. 

To determine their local density, not only must wimps be observed both in direct detection experiments 
and at colliders, but it must be established that these two signals correspond to the same particle. The only 
known way to suggest that the same particle is observed in both detectors and colliders is to compare the 
masses of the particles seen in each. This requires that the wimp mass be determined independently by both 
direct detection experiments and colliders. We reviewed the heretofore known way to determine the wimp 
mass using direct detection experiments, and presented a somewhat new, explicit technique. We discussed 
how, in contrast to direct detection experiments, there does not exist any general method of determining 
the mass of wimps produced at hadron colliders, though it may be possible in particular lucky cases. 

To highlight the general arguments presented, we often referred to the case of supersymmetric dark matter, 
where the discovered wimp is the LSP. We criticized the trend of studying dark matter by calculating the 



If the gauge content of the neutralino was unknown, the interaction parameters could have been maximized with respect to 
these parameters as described earlier. Therefore, the upper bounds obtained are more restrictive than would be in practice. 
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thermal freeze-out contribution to the LSP relic density. We showed that determining even the (quintessence- 
free) thermally-produced contribution to the relic density of the LSP is both fundamentally insufficient and 
likely a practical impossibility using data from hadron colliders. In general, too many (possibly) difficult to 
measure lagrangian parameters are required to compute this for a generic supersymmetric model; and even 
if all the parameters were known, the result can be extremely sensitive to high-scale physics unobservable at 
colliders. 

In contrast to thermal freeze-out calculations, we found that the LSP-nucleon incoherent scattering cross 
section could be absolutely bounded (and perhaps well-approximated) using only sparse data from colliders 
for the most general MSSM (without any presumptions about supersymmetry breaking); furthermore, we 
showed that this bound is iteratively improvable as more data becomes available from colliders. This bound 
was shown to be robust for thousands of realistic test models. Combining this bound on the scattering cross 
section with data from direct detection experiments immediately allows one to obtain into an absolute lower 
bound on the local density of neutralino dark matter. 

This is the first study of which we are aware that fully addresses the difficulty of determining what fraction 
of dark matter is accounted for by wimps discovered at colliders or in detectors. We have shown explicitly 
how to deal with this problem in the case of supersymmetric dark matter, and described how similar analyses 
could be done for any other dark matter particle candidate. Although the dark matter problem may not 
be solved immediately when wimps are first observed, there are clear and general ways to address their 
cosmological significance. 
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